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WEAK DUAL PAIRS AND JETLET METHODS FOR IDEAL 
INCOMPRESSIBLE FLUID MODELS IN n > 2 DIMENSIONS 

C. J. COTTER 1 , J. ELDERING 1 , D. D. HOLM 1 , H. O. JACOBS 1 ’*, AND D. M. MEIER 2 


Abstract. We review the role of dual pairs in mechanics and use them to derive particle- 
likc solutions to regularized incompressible fluid systems. In our case we have a dual pair 
resulting from the action of diffeomorphisms on point particles (essentially by moving the 
points). We then augment our dual pair by considering the action of diffeomorphisms on 
Taylor series, also known as jets. The augmented weak dual pairs induce a hierarchy of 
particle-like solutions and conservation laws with particles carrying a copy of a jet group. 
We call these augmented particles jetlets. The jet groups serve as finite-dimensional models 
of the diffeomorphism group itself, and so the jetlet particles serve as a finite-dimensional 
model of the self-similarity exhibited by ideal incompressible fluids. The conservation law 
associated to jetlet solutions is shown to be a shadow of Kelvin’s circulation theorem. Finally, 
we study the dynamics of infinite time particle mergers. We prove that two merging particles 
at the zeroth level in the hierarchy yield dynamics which asymptotically approach that of 
a single particle in the first level in the hierarchy. This merging behavior is then verified 
numerically as well as the exchange of angular momentum which must occur during a near 
collision of two particles. The resulting particle-like solutions suggest a new class of meshless 
methods which work in dimensions n > 2 and which exhibit a shadow of Kelvin’s circulation 
theorem. More broadly, this provides one of the first finite-dimensional models of self¬ 
similarity in ideal fluids. 


1. Introduction 

Arnold’s geometric insight in |Arn66| has forever changed the way mathematicians look 
at ideal fluid dynamics. According to |Arn66| . ideal incompressible fluid motion on an 
orientable Riemannian manifold M is equivalent to geodesic motion on the Lie group of 
volume preserving diffeomorphisms SDiff(M) (i.e., the smooth invertible volume preserving 
maps of M into itself, with smooth inverses). The Riemannian metric is simply the fluid 
kinetic energy, which is the L 2 -norm of the fluid’s velocity field. This characterization of 
ideal fluid flow allowed Poisson geometers and geometric mechanicians to provide a new 
perspective on ideal fluids and other PDEs with hydrodynamics background [EM701 IMW83I. 
IZeiQll IAK981 IHMR981 IFHTOlj . As will be shown, such a perspective is particularly fruitful 
for the purpose of reducing the infinite-dimensional fluid system to a finite-dimensional 
ordinary differential equation (e.g. point vortex solutions). 

The dimension reduction that will be performed in this paper and the method which it 
proposes bears much semblance to the point vortex method [ Cho73] . In particular, ' MU S3 
illustrated how the point vortex solutions and the conservation of circulation could be derived 
via a pair of Poisson maps known as a dual pair. One of these Poisson maps was the 
embedding map from a finite-dimensional manifold into the infinite-dimensional space dual 
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to the divergence free vector fields. In other words, the space of point vortices is merely 
an invariant manifold of ideal fluid motion]]] It was later realized that the use of dual pairs 
related to fluid applications was problematic and needed to be relaxed. This led to the 
notion of weak dual pairs [ GBV12) . In this paper we will derive a hierarchy of different weak 
dual pairs, in order to obtain a class of hnite-dimensionally parametrized solutions, each of 
which comes with a conservation law that shadows Kelvin’s circulation theorem. Just as one 
can consider the point vortices to be the atoms of the point vortex method, the atoms of 
these new solutions are particle-like objects which we call jetlets (or k-jetlets if we wish to 
emphasize that we are considering the k -th level of the hierarchy). 

It is notable that the zeroth level of the hierarchy is a classical particle-like solution which 
appears in many geodesic systems on diffeomorphism groups ( CH93 , JM001IFH011 lMM13bj . 
It is also notable that, unlike a point vortex, a jetlet is also well-defined in dimensions greater 
than two. 


1.1. Main contributions. In this article we derive a hierarchy of particle-like solutions for 
a regularized model of ideal fluids described in [MM 13b) as motion on the group of volume 
preserving diffeomorphisms of M n , denoted SDiff(R n ). Each level in the hierarchy consists 
of particles with internal group variables, parametrized by a finite-dimensional model of a 
diffeomorphism group known as a jet group [KMS991 Chapter 4], The jet groups lie at 
the foundations of certain representation theories of diffeomorphism groups (see [ VGG75I 
Appendix 2] or | lKir81j ). Hence, it seems natural to invoke these foundations in the context 
of fluids, where the configuration manifold is a diffeomorphism group. The particle-like 
solutions, which “carry” jet groups, paint an intuitive picture of a large scale diffeomorphism 
that advects particles, each of which carries its own partial description of a local deformation 
of the fluid in a small region around it. At each higher level in the hierarchy, the description of 
the deformation becomes more detailed. Thus the jets of diffeomorphisms possess a natural 
sense of the “self-similarity” which is present in the diffeomorphism group itself. Models 
such as this are crucial to our understanding of fluids, both from a numerical perspective 
and from the perspective of fundamental mathematics. 

Specifically, we will accomplish the following: 

(1) We provide the first explicit Hamiltonian description of the full jet hierarchy of the 
particle-like solutions discovered in [ DJR13!. GHJM14] . 

(2) We compute a nested sequence of conserved quantities at each level in the hierarchy. 
Each of these conserved quantities will be related to the conservation of circulation. 

(3) We numerically compute some particle-like solutions at the zeroth and first levels in 
the hierarchy. We will observe a form of cascade in which interactions of solutions at 
the k-th level tend asymptotically in time toward solutions at the (k + l)-th level. 


1.2. Outline of the paper. After the introduction, in Section [2j we outline our strategy: 
The Erst goal is to show that jetlet models admit a weak dual pair at each level in the 
hierarchy. Once this has been done, Theorem A.9 (proven in the appendix) does the rest. 
Namely, the jetlets satisfy canonical Hamiltonian equations, and the right momentum map 
is conserved by the flow. 


1 This is not quite correct, as one must ignore the infinite self-energy terms of point vortices and then 
extend the space of admissible solutions to non-smooth velocity helds. However, modulo this physically 
motivated caveat the statement holds. 
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Section [3] contains a brief discussion of the standard dual pair for the example of the rigid 
body. 

In Section [4] the main development of the paper starts. We first discuss briefly the Lie- 
Poisson approach to ideal fluids and Euler’s equation. This is to set the scene for what 
follows. In Section 4.1 we discuss the Mumford-Michor model |MM13b] and mention its 
standard dual pair (in parallel with the rigid body, the two legs of the dual pair correspond 
to the cotangent lift momentum maps for right and left actions, respectively). We also 
recall that the conservation of the right momentum map Jr here is equivalent to Kelvin’s 
circulation theorem. Again, the dual pair viewpoint is extremely efficient. Once one has 
realized that there is a dual pair, one knows (from right-invariance) that the left momentum 
map Jr maps Hamilton’s equations on T* SDiff to (reduced) Lie-Poisson form, and at the 
same time that Jr is conserved. 

In Section 4.2 we introduce zeroth order jetlets, also called landmarks, and we discuss 
their dual pair. This is done in parallel with later sections. The left momentum map is the 
usual one, and the right momentum map is in fact trivial. So the ‘dual pair’ is somewhat 
unnatural here, but it helps present already the kind of thinking that will be employed in 
the later parts of the paper. 

In Section 43 we consider Erst order jetlets. Everything here goes in parallel, except that 
one has to be more careful when introducing the relevant right and left actions. Once one 
has the definitions in place, it is not difficult to recognize that the cotangent lifts for the right 
and left actions lead to, now, a weak dual pair of momentum maps. We have to add ‘weak’ 
here since the transitivity of the left action on level sets of the right momentum map is lost, 
but a weak dual pair is retained, essentially because the group actions still commute, see 
Definition A.8 Proposition |4.5| collects the results that follow immediately as a consequence 
of the weak dual pair. 

Throughout the paper we have taken care to explain the intuition behind the more abstract 
concepts. The discussion just after Proposition 45 is an example of this. There, we discuss 
the relationship between the jetlet solutions on the diffeomorphism group and the resulting 
trajectory on the space of Taylor jets. 

In Section 4.4 we describe jetlets at general levels. Again, once the relevant spaces and 
actions have been introduced (which is now quite an intricate endeavor), one recognizes that 
the momentum maps are a weak dual pair by the usual arguments. We give the form of 
the momentum maps explicitly just after Proposition |4.6[ Again, the weak dual pair leads 
to analogous conclusions about dynamics. Namely, it is canonically Hamiltonian, with Jr 
conserved; Jl maps the dynamics to the Lie-Poisson dynamics on the one-form densities. 

In Section 45 we discuss in more detail Kelvin’s circulation theorem. More precisely, 
we discuss the relationship between the ‘standard’ circulation theorem for the fluid and the 
conserved momentum maps Jr at the various levels of jetlets. The main result is represented 
schematically and proven in detail (all that is required, in essence, is the standard formula for 
cotangent lift momentum maps). Intuitively speaking, what we show is that the conserved 
jetlet momentum maps are ‘shadows’ of the ‘full’ right momentum map of the fluid. 

In Section [5] we discuss the explicit dynamical behavior of the particle model, in particular 
we study ‘collisions’ leading to mergers of jetlets. This makes an explicit connection between 
dynamics in the different levels of the jetlet hierarchy. It can potentially be useful when 
simulating jetlet systems: when two particles become close, they can be replaced by a merged 
state whose momenta do not blow up. 
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We redo the analysis of Mumford and Michor to find that two 0-jetlets can either merge 
in infinite time, or bounce off each other. Then, we analyze the asymptotic dynamics of this 
merged state and show that it coincides with the dynamics of a single 1-jet let particle. This 
improves the claim in (CHJM14 J which showed the convergence to a 1-jetlet state without 
explicitly considering the dynamics. 

Finally, we suggest a more algebraic interpretation of mergers in the jetlet hierarchy by 
viewing all levels of the hierarchy as embedded in the larger space Xdi v (18C)*, where levels 
form (part of) the boundaries of other levels (a bit like in a CW-complex). 

In summary, the numerics section [6] shows first of all that it is feasible to implement this 
jetlet model numerically, also in dimensions higher than two. Further, we corroborate the 
analytical results and confirm the conserved quantities and jetlet particle merging behavior 
numerically. Several different experiments show that the merging/scattering behavior per¬ 
sists under various perturbations that cannot be studied analytically anymore. It also shows 
how Jl (i.e. angular and linear) momentum is exchanged in the jet-particle collisions. 

We provide a detailed appendix. In Appendix [A] a brief discussion of symplectic and 
Poisson manifolds is followed by the definition of weak dual pairs. Then we prove Theo- 
|X9l 


rein 


which is the core result used in the paper. Appendix [B] presents 


an overview over 


the spaces used in the paper and relates them to some general results of reduction theory, 
see in particular Figures [6] and [7J In Appendix [C] we describe the index conventions used in 
the main text (e.g., when calculating momentum maps for the general jetlet solutions), while 
Appendix [D] briefly discusses the more abstract point of how the dual space of vector fields 
can be viewed as tensor product of 1-forms and distributions. Finally, Appendix [E] provides 
information describing the equations of motion for 1-jets, but in reduced coordinates, which 
arise after the reduction that eliminates the conserved quantities, Jr. 


1.3. Previous work. Lagrangian models of ideal fluids such as smooth particle hydrody¬ 
namics (GM77 . ILuc77 j and vortex methods [ Cho731| . do not exhibit structures which ex¬ 
press the nested character of the diffeomorphism group. One means of obtaining a La¬ 
grangian model with a nested structure was recently presented in [DJR13j where a sequence 
of infinite-dimensional reductions by symmetry was executed, to produce a hierarchy of 
finite-dimensional systems from a regularized fluid model. The finite-dimensional systems 
were particle-like solutions in which each particle carries a model of the diffeomorphism 
group, known as jet groups. Thus, [ D.TR13 ] derived a Lagrangian analog of “whirls within 
whirls”, similar to the Eulerian models proposed in jHT12j . However, the specifics of the 
regularized fluid model were not determined, and the analysis of DJR13 - was purely formal. 

Later, a regularized version of the ideal fluid equations was presented by (MM13bj . This 
new partial differential equation was amenable to the theory presented in I D.IR 13 . and 
gave rise to a specific and easily implementable manifestation of the hierarchy of particle 
models described there. In )CH.TM14 ] we numerically computed some of these particle-like 
solutions and observed cascade phenomena as an emergent behavior at the zeroth level in 
the hierarchy. 

It is notable that the zeroth level of the hierarchy has been studied in the context of partial 
differential equations with hydrodynamic background. In particular, [ HR06j provide the first 
convergence proof of such a method in the context of the Camassa-Holm equation. This same 
idea was implemented for the n-dimensional Camassa-Holm equation in 1CDTM12] . In the 
context of image registration algorithms, the need to obtain compressible diffeomorphisms 
motivated the use of particle methods in a similar framework jJMOO j. These methods, 
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designed for a wide array of PDEs, are studied analytically in [ TY05j . which also contains 
a proof of wcll-posedness for a range of PDEs. Finally, in the context of image registration, 
[ SNDP13 ] discovered a compressible fluid version of the hierarchy derived in .D JR13 I, and 
numerically integrated solutions in the first level of the hierarchy. 

1.4. Notation. We will let £(M n ) denote the space of H°° vector fields, and we let Xdi v (BP) 
denote the space of divergence free vector fields resulting from the Hodge decomposition. 
We let Diff(M n ) denote the space of H°° diffeomorphisms of M n (see [MM13aj ) and we let 
SDiff(M n ) denote the subgroup of volume preserving diffeomorphisms. 

Various different types of indices will be used throughout this paper. To distinguish 
between the types, we keep the following conventions: 

• indices a, 6 , c,... label particles and range from 1 to N] 

• indices i,j, k,... label space coordinates and range from 1 to n, the dimension of 
space; 

• superscript indices ( k ) denote the order of jets in a coordinate-free representation. 


2. Main approach 


Let us begin by describing our strategy for obtaining particle-like solutions to regularized 
fluid equations. Our main framework is that of Hamiltonian mechanics and symplectic 
geometry. In particular, our main hammer is Theorem A.9 (see page 31 in the Appendix), 
which we repeat here for convenience. 


Theorem A. 9 Let Pi and P ' 2 be Poisson manifolds and let S' be a symplectic 
manifold. Let Ji, J 2 : S — » Pi, P 2 form a weak dual pair, see Definition |A.8 
Let h G C 1 (Pi). If (g,p)(t) G S is a solution to Hamilton’s equations with 
respect to the Hamiltonian H = h o J 1; then Ji ((g,p)(t)) G Pi is a solution 
to Hamilton’s equations on Pi with respect to h, and ((<?,£>)(£)) is constant 
in time. 


We will leverage this theorem in the following way. We will find a (weak) dual pair, 
Ji, J 2 : S —> 0 *, where 

(1) S is the space of particle locations and momenta (and, for the higher orders of the 
hierarchy, internal group variables), 

( 2 ) g is the space of vector fields, 

(3) and h is the Hamiltonian of a regularized model of ideal fluid (i.e. a kinetic energy). 

Theorem [A9] applied to the (weak) dual pair Ji, J 2 then yields particle-like solutions and con¬ 
served quantities to the equations of motion of a regularized model of an ideal incompressible 
fluid. 


3. Example: the rigid body 

In this section we review basic notions from classical mechanics by studying the motion 
of a rigid body whose center of mass rests at the origin. We will see a first application of 
Theorem A.9 in this context. The configuration of a rigid body is described by a rotation 
matrix R G SO(3). The equations of motion are given by Hamilton’s equations on the 
cotangent bundle T*SO(3). There are canonical coordinates on T* SO(3) given by (P, P) 
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where P is such that R T P is a 3 x 3 anti-symmetric matrix^ The angular momentum in the 
body frame is the unique vector n G R 2 3 = so(3)* such that 


Jr(R,P ) : =R T P = 


0 —II 3 II 2 

n 3 0 “iii 

—1I2 iii 0 


Denoting this map from (R, P) to II by Jr, the Hamiltonian can be written as a function of 
n. In particular, the reduced Hamiltonian is 

K n) = ^n-r 1 -n, 

where I is a non-degenerate 3x3 symmetric matrix known as the moment of inertia matrix. 
The unreduced Hamiltonian is H(R,P) = h(J r(R, P)). 

The group SO(3) acts upon itself by left multiplication. This action can be lifted to an 
action on T* SO(3) given by 


(R, P ) G T* SO(3) y TT ' {gR, gP ) G T* SO(3). 

It is notable that n = R T P = R T g T gP = ( gR) T (gP ) is unaltered by this transformation. 
Therefore the Hamiltonian, /i(n), is invariant under this action of SO(3). By Noether’s 
theorem there is a conserved quantity associated to this symmetry. The conserved quantity 
is manifested by the momentum map 

J l (R,P) = PR t . 


Moreover, due to this symmetry, one can write the evolution equations on a lower-dimensional 
space. The very fact that the Hamiltonian is written in terms of n suggests that the equations 
of motion can be written in terms of n alone. Indeed this is the case, 

(1) ri = n x (i^n). 


This equation can be seen as a Hamiltonian equation on R 3 with respect to the non-canonical 
Poisson bracket 


{F, G} Nambu (t) = —X -(VFX VG) 


known as the Nambu bracket. This is no coincidence. Let us first introduce the isomorphism 
of Lie algebras (R 3 , x) and so (3), given by the so-called hat map 


( 2 ) 


x G R 3 1 —>■ x 


0 -x 3 x 2 

x 3 0 —x\ 

—x 2 Xi 0 


G so(3). 


Proposition 3.1 (see [Hollli §2.5]). The Nambu bracket on R 3 is identified with the Lie - 
Poisson bracket onso(3)* through the hat map isomorphism (|2]) in the sense that {/, gjNambu^ 
{f,9} lp(£) where f,g G C^R 3 ) and f,g G C' 1 (so(3)*) are defined by f(x) = f(x),g(x) = 
g(x). 


2 We use the pairing (P,v) = tr (P T v) for P £ T^SO( 3) and v £ TrSO( 3). More generally, the pairing 

(■, ■) will always refer to the natural pairing between a vector space and its dual in this paper. 
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Proof. The Lie-Poisson bracket on so(3)* is 


{f,gh p(n) 


n. 




for arbitrary functions f,gE C' 1 (so(3)*). There exist functions f,gE C 1 (M 3 ) related to F, G 
through the hat map. One can observe that df(fl) E s o(3) is related to V/(fI) through the 
relation V/(L1) = d/(fl). We see that the commutator bracket satisfies 


[x, y\ := xy — yx = x x y. 


Therefore the Lie-Poisson bracket can be written as 

= (n, v/ITv^n)) 
= n ■ (v/ X Vg) (n) 

{/;!?} Nambu (II). 


□ 


Recall that T* SO(3) is a symplectic manifold. The momentum maps Jl, Jr arise canon¬ 
ically from the left and right action of SO(3) on T*SO(3). The actions commute and it 
can also be checked that Jl and Jr have symplectically orthogonal kernels, hence it follows 
from jGBV12l Corollary 2.6] that the diagram 


so(3)* -f-k- T* SO(3) ^%so(3)* 


PR 


T Jl 


(R.P) 


Jr, 


r 2 p 


is a dual pair (see page 31 for details). The maps Jl and Jr are called symplectic variables 

S as they allow one to pull-back calculations on a Poisson manifold to a symplectic 


m 

manifold 


By Theorem A.9 this dual pair expresses rigid body dynamics and conserved quantities. 


Specifically, the right leg yields the reduced phase space where the system evolves in time. 
The left leg yields the conserved quantities of the rigid body associated with the left action 
of SO(3) on itself. The most important aspect of these maps is that they are both Poisson 
maps, i.e. they carry the canonical Poisson bracket on T* SO(3) to the Nambu bracket on 
so(3)* = (M 3 , x), as the following proposition shows. 


Proposition 3.2 (remark 2.5.11 (Holllj ). Let {•, -} cari denote the canonical Poisson bracket 
on the cotangent bundle T* SO(3). Let {•, -}Nambu denote the Nambu bracket on M 3 . Both J L 
and Jr are Poisson maps. Explicitly, this means 


{f ° Jli g ° ^L}can — {/; S^Nambu 0 Jl 

{/ ° Jri g ° Jr \can {/> g }Nambu ° JR 


for any f,g E C 2 (M 3 ). 


3 , MW83( also referred to Jr and Jl as Clebsch variables, however this terminology has changed over the 
past few decades. 
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Proof. Let (R, P ) G T* SO(3) and set II = Jr(R, P). We observe that 

(n, ft) = tr(n T ft) = tr (P t rCi) = (P,R-Cl). 

This tells us that Jr is the momentum map associated with the cotangent lift of the right 
action of SO(3) on itself. Such momentum maps are always equivariant and thus yield 
Poisson maps (Theorem 12.4.1 [MR99]). Thus Jr carries the canonical Poisson bracket on 
T* SO(3), to the Lie-Poisson bracket on so(3)*. By Proposition 3.1, this is nothing but the 
Nambu bracket upon identifying so (3) with M 3 . The same argument applies to Jl using a 
left action. □ 

One can obtain solutions to ([TJ) by solving canonical Hamiltonian equations with respect 
to H(R,P). In particular, if (. R,P)(t ) is a solution to Hamilton’s equation, then n(t) = 
Jr((R , P)(t)) is a solution to Hamilton’s equation with respect to the Nambu bracket. This 
is a result of Proposition A.7| paired with the observation that Jr is a Poisson map via 
Proposition |3.2 


4. Regularized fluids 


Euler’s equations of motion for incompressible fluids can be seen as Hamiltonian equations 
on the (dual) space of divergence free vector fields [ ArnGfij . Consider the Lie algebra of vector 
fields on M n , denoted by X(M”). Formally, the dual space to j£(M n ) is a Poisson manifold 
when equipped with the Lie-Poisson bracket (see (27) in Appendix A.2). We may consider 
the map 0 : T*R n —* £(M n )* given implicitly by 

(0(g,p), u) = p ■ u(q ) Vw G £(M n ), (q,p) G T*R n . 

Explicitly we may write 0 using the Dirac-delta functional as 0(g,p) = p®S q . It is shown 
in [HMQ5J that this map is Poisson. Furthermore, as the n-dimensional Camassa-Holm 
equation [ICH93] is a Hamiltonian equation on j£(M n )*, Proposition A.7 promises to express 


a certain subset of solutions by solving Hamiltonian equations for a finite number of particles. 
Specifically, 0 yields the peakon solutions of the n-dimensional Camassa-Holm equation. In 
this section we explore analogous constructions for an incompressible and regularized version 
of the Camassa-Holm equation, discovered in |MM13bj . 

In the case where /iRuier(rn) = |||m ||^ 2 is the standard fluid kinetic energy on the dual 
vector space to the incompressible vector fields, SdivCR")*, Hamilton’s equations are written 


as 


(3) d t m + C u [rn] =0 , u l = 8 l ^mj. 

where C u [m\ is the Lie derivative of m. The primary finding of (Arn66j was that (|3]) is 
equivalent to the inviscid fluid equation 

d t u + u ■ Wu = — Vp , V • u — 0. 

Since we have not yet clarified the Poisson structures of the system, it may not be obvious 
that (J3]) is a Hamiltonian equation. The following proposition shows this for a general 
Hamiltonian on Xd; v (M n )*. 

Proposition 4.1 (I Ani66 ). Let h G C°°(3Cdiv(BC)*)- Recall that 3C d iv(BC)* Is a Poisson 
manifold when equipped with the Lie-Poisson bracket, and given a function h, the Frechet 
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derivative dh(m ) is an element o/£ d i v (M n )**. In the event that dh(m ) G £ d i V (M n ), Hamilton's 
equations are given by 

fn + £ u [m] =0 , u — dh(m). 


Proof. To each v G £ d i V (M n ) we can associate a linear function on £ div (M n )* given by m G 
£ div (M n )* i-> ( m,v) G M. Let us denote this function by /,,. Let m{t) satisfy Hamilton’s 
equations. By the dehnition of the Lie-Poisson bracket (see (27) in Appendix |A.2[ ) we observe 

= {fv,h}{m) = (m, [df v (m),dh(m )]) 

= (m,jCdh(m)[dfv(m)}). 


However df v = v and so the last line can be equated with — (£ dh ^[m],v). Additionally, we 
know that 4 f v (m) = {rh,v) since v is constant in time. Therefore we find 

(m + £ d h(m)[m],v) = 0 . 

As v is arbitrary, this uniquely characterizes ?h|^] The result follows. □ 


4.1. The Mumford—Michor model. Consider the Hamiltonian 

h P} a(m) = ^ (m, K P:(T * m) L 2, 

where K Pt<T : M n —> M. nxn is the matrix valued Green’s kernel defined by the property 

(l - —■ / K^Jx - y)m j {y)dy = d^m^x). 

V P y 

In this case Hamilton’s equations take the form 

(4) d t m + £ u [m\ =0 , u l = K^ a * rn r 

Solutions to Q exhibit existence and uniqueness for all time. Moreover, as a —> 0, h pa -A 
hEuler and one can speculate that solutions to Q approach solutions to the ideal fluid equa¬ 
tion (g. In fact, this is the case over short times, and for a > 0 solutions of Q differ 
from those of (|3| by an amount at in the H k - norm. Thus Hamilton’s equations with respect 
to h P:Cr have been proposed as a model for ideal fluids [ MM13b i Theorems 2 and 3]. From 
now on we shall often suppress the parameters p, a and shorten K = K p a to prevent index 
clutter. 

Next, we discuss a dual pair for this system. The group SDiff'(M n ) acts on itself from the 
left and from the right. These actions can be lifted to T* SDiff(M n ), and yield momentum 
maps Jl , Jr - T* SDiff(M n ) —s- £ div (M n )*. In particular these maps form the dual pair 

£ div (M n )* g-Lh T* SDiff(M n ) 4% £ div (M n )L 


For a mathematically rigorous treatment of this dual pair we refer to |GBV12j . 

By Theorem A.9 we can use this dual pair to derive dynamical properties of Hamiltonian 
equations defined on £ d i v (M n )*. The Hamiltonian for a fluid is written on the left instance of 
X d ; v (M n )*. One can (in principle) solve Hamilton’s equations on T* SDiff(M n ) with respect to 
the Hamiltonian H = h o J L . This yields the material or Lagrangian coordinate perspective 
of fluid mechanics one encounters in a first course on continuum mechanics. The right leg 


4 This is not a “weak” characterization. The entity to is contained in the dual space to Xdi v (R n ) and it is 
therefore defined uniquely by how it acts on Xdi v (K"). 
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yields conserved quantities associated with the particle relabeling symmetry of the fluid. It 
was found in |Arn 66 ] that these conserved momenta are identical to the law of conservation 
of circulation, that is, Kelvin’s circulation theorem. 

Unfortunately, this dual pair does not help us in solving Q since solving Hamilton’s 
equations on T* SDiff(M n ) is no less difficult than solving Hamilton’s equations on Xdi V (®U)*- 
In the next section we will derive a dual pair wherein the symplectic manifold is more 
reasonable. This will yield the particle-like solutions described in [IMM13b| . In the later 
parts of the paper we will generalize the treatment to obtain a weak dual pair for each level 
in the hierarchy of particle-like solutions. 


4.2. Particle-like solutions. There is a natural left group and algebra action of SDiff(lk 
and 3Cdiv(K n ), respectively, on M n given by 

^eSDifRir) 


q e 


ue.L'di, 


<p(q) G IP 
’ u{q ) G TqR n . 


? 6 ^ 1 - t uynj sc J-q 

The tangent lift of the former is defined in the obvious way by sending 

(q,v) G TW 1 = K n x I” (d^iq^di) G T v{q) R n . 

The cotangent lift is defined by taking the dual of the tangent lifted action. That is to say, 

¥ jeSDiff(R , q 


(q,p) e T* IP = IP X 


(<p 1 (q),(d i ip 3 (q)p j )e l ), 


where {e*} forms the dual basis to {<9*} at p 1 (g). The momentum map J : T* 
Xdiv(K n )* associated to this left action is defined by the condition 

(Jl\<LP),u) ■= (pX<?)), 


-A 


for all u G Xdi 


divl 


*) and (q,p) G T*M. n (see (26) in Appendix A.l); the superscript in J 


(o) 


serves as a reminder that we are considering the zeroth level in the hierarchy of particle-like 
solutions. We see that J^\q,p) is an evaluation operator, and we can write it more explicitly 
as a measure-valued momentum map in terms of the Dirac-delta distribution as 

Jl\^p) =P®$qi 

where S q is the Dirac-delta distribution on M n centered at q. This identification holds modulo 
dP 1 (M n ) ig) dx where dC 1 (M r) ) is the space of C°-exact one-forms and dx is the canonical 
volume form on (see Appendix D). Since is a cotangent lift momentum map, it is 
equivariant and therefore Poisson, again by ( MR99 1 Theorem 12.4.9]. 

We define the manifold for particles, 

Q‘k 


'N ~ {(?!) ■ ■ ■ ) Qn) £ 


x • • • x IP | q a 7 ^ q b when a ^ b} . 


We will index the particles with a, b, c,... and Cartesian coordinate directions in M n with 
indices i,j, k,... Thus each q e can be decomposed into N particles as (q x ,..., q^), 
where each q a G M n and the f-th coordinate of the a-th particle is denoted by q a \ 

The group SDiff(M n ) acts on by the diagonal action. Through the same manipulations 
as we applied previously we obtain the momentum map for N particles given by 

4 0) (? ; P) =Pa®fiq a , 

where a sum over repeated indices is implied. 
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Figure 1. A 0-jetlet with momentum m = (1, 0) 7 ® d 0 in dimensions 2 and 3. 
(left) Quiver plot of the induced 2 dimensional velocity field, (right) Streamline 
plot of the induced 3 dimensional velocity field. 


By Propositions A. 7 and 4.1 we obtain solutions to Hamilton’s equations on j£ div 
by solving Hamilton’s equations on T*M n if dh(J l f\q, p)) is a vector held 
we calculate that dh evaluated on the J^ (p, q) is the vector held 


If h = h p ai then 


dh(p a ® S 9 a )(x) = K lJ (x - q a )Paj-Q-- 

This is a vector held whose differentiability is determined completely by that of the ker¬ 
nel, K. Once one has found a solution (q(t),p(t)) to Hamilton’s equations on T*M. n and 
thus also a solution to Hamilton’s equations on £ div (M n )*, one can proceed to integrate the 
corresponding time-dependent vector held to obtain the fluid motion <p t in SDiff(M n ). It is 
natural (but not mandatory) to choose the initial map <po to be the identity or at least to 
be a diffeomorphism that satishes <£o(<7a(0)) = q a ( 0) for all a. If this choice is made, one can 
interpret the curve q(t) = (qi(t),... , gv(Q) in Qn as the locations of particles as they are 
swept along by the fluid how, that is, q a (t) = Pt(q a ( 0)) for all a — 1,..., N. 

We see that J)° is injective, and thus has a trivial kernel. As a result, the symplectic 
orthogonal to the kernel of is the full tangent bundle T(T*Q$). Hence, if we define the 
(trivial) map : T*Q$ —> 3C div (M n )* by = 0 , it follows that the diagram 


^(MT 

AO) 

•fL 

T*Q ( n ] 

r(0) 

J R ^ 

X div (M n )* 


j' 0) 


J (° } 


Pa ® Sq a 

^ L j 

(q,p) 

1 .R ^ 

0 


is a (proper) dual pair. This dual pair allows us to express conservation laws and dynamics 
as a result of Theorem |A.9| Namely, the left leg represents the space in which particle-like 


(o) 


solutions to Q evolve, while the right leg represents a (trivial) conserved quantity. 

In order to make contact with the later parts of the paper, it is useful to remark that 
can formally be understood as the cotangent lift momentum map associated with a certain 
(trivial) group action. To that end, we fix a designated point z = (zi, ..., Z/y) G and 
take an arbitrary element q = (qi, ... , qjy) G to represent the set of all ip G SDiff(M”) 
that satisfy <p ■ z = q. That is, p(z a ) = q a for all a. This means in particular that the 
specification of (qi,..., q n) fixes the zeroth order Taylor expansion (at the locations z a ) of 
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the corresponding set of diffeomorphisms. With this in mind, let us define the isotropy group 

(5) iso(^) = {if G SDiff(M n )|-0(^ a ) = z a for all a} 

and a (trivial) right action on where if £ iso(z) maps an element (p(zi),..., p(zn)) 
to (p o 'ijj(zi ),..., ^ o if(zjsr)). Similar constructions will be crucial in later sections, when 
constructing the right leg of an — in this case weak — dual pair for the higher levels in the 
hierarchy of particle-like solutions. 

In summary we find: 

Proposition 4.2 (§7 of | IMM13bj ). Let p > | + 1 and a > 0. Let H : T*Q$ — > R be the 
function 

H(q,p ) = ^KpAq a ~ qb)PaiPbj- 

If (q,p)(t) is a solution to Hamilton’s equations, then m(t ) = p a {t) ® 5 q ( t ) is a solution to 

0 - 


for more details on the kernel smoothness condition on p. □ 


Proof. 


Note that H = h Pi(T o J L 


( 0 ) 


and apply Theorem 


A.9 


See the remark after Theorem 


4.8 


4.3. First order particle-like solutions. In this section we revisit the first order particle¬ 
like solutions of [CHJM14] and discuss their weak dual pair, before extending the treatment 
to the higher levels of the hierarchy in the subsequent section. Let SL(n) denote the Lie 
group of n x n matrices with unit determinant. Let q = (q^°\ q dl) G R n x SL(n) and consider 
the left SDiff(R n ) action on R n x SL(n) given by 

(6) p ■ q = (p(q {0) ),Dp(q (0) ) ■ q W ). 

Where Dp(q^) ■ gd) j s the result of multiplying the Jacobian matrix Dp(q^) with gd). 

Proposition 4.3. The action o/SDiff(M) on R x SL(n) in @ is a group action. 

Proof. Since p G SDiff(R n ), it follows that Dp | g (o) G SL(n). Therefore Dp\ q (o) • gd) G SL(n). 
Secondly, if pi,p 2 £ SDiff(R n ) we observe that 

P2 ■ (<Pi -q) = (p 2 • (<Mg (0) ), D<Pi\ q (o) ■ q (1) ) 

= D T 2 \ Mq ( 0)) ■ Dp i| g( 0 ) • g (1) ) 

= ((a o a)(9 (0) ).%2°7i) 1,(0) -q (1) ) 

= (P 2 °pi) -q, 

where the second and third lines are applications of the chain rule. □ 


As before, this action can be lifted to the cotangent bundle T*(R n x SL(n)). Specifically, 
the action of the diffeomorphism p -1 is given by 

(7) (7»* • (<I {0) ,q m ,P l0} ,p m ) = [Dp I ,,.,]" 1 ■ q m ,Dp\' qa) ■ p m , Dp |J„, 

Also as before, we can generalize this construction to the space of iV > 1 particles by 
considering the space 


Q ( n = (Qu---,Qn) 


= {qa\qa ] 

g(°) -/ AO) 




) gK"x SL(n) 
when a ^ i 
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For convenience it is nice to choose coordinates at this point. If we let q a denote the position 
of the a-th particle, q a l denote the i-th component of this position, and let q a t J denote the 
(i,j) entry of the a-th matrix then the resulting momentum map is defined by the condition 

(Jl\q,p),u) =p ai u\q^ ) )+p ai 3 d k u\q ( £ ) )q a k j 

for an arbitrary u G 2 div(M n ). In terms of the Dirac-delta functional, we can write as 
the measure-valued momentum map 

Jl\q,p) = PaA xl ® 8 q ( 0 ) - pjq a k jdx' ® d k S q ( 0 ). 

Next, we construct a dual momentum map associated with a right action on Q$. Anal¬ 
ogous to the previous section it is useful at this stage to write elements of Q$ in the 
form p ■ z = (p(z\, 1),..., p ■ (zjv, 1)), p G SDiff(M n ), for some designated element z = 
((zi, 1),..., (zjv, 1))- Clearly, every element of can be written in this form for some 
p G SDiff(M"). Indeed, with this convention the specification of an element in hxes the 
first order Taylor expansion of p at the locations z a , a = 1,..., N. With this in mind, let us 
recall the isotropy group iso(z) defined earlier in (|5]) , which leaves these locations invariant, 
and define a right action on given by 

(p ■ z) ■ if — (p o ip) ■ z. 


That is, 

( 8 ) 


( 9 < 0 ) , 9 (1) 


)A=(( 


qi\qi ] 


DiP\ 


z 1 


( 0 ) ( 1 ) 
Qn 


Di>\ 


Z N 


In Proposition 4.6 we will generalize this construction to define a right action for the higher 
levels in the hierarchy of particle-like solutions. The action (J 8 | yields the cotangent lift 
momentum map defined by the condition 

(Jr\q,P),u) = Pai(lak d j uk ( Z a)- 


In terms of the Dirac-delta functional we may write this as the measure-valued momentum 
map 


(9) 

Proposition 4.4. 


Jr VhP) = ~PaiQak ( 
The momentum maps J^ 1 and 


form a weak dual pair. 


We postpone the proof, as this is a special case of a proposition which comes later in the 
paper (Proposition |4.7[ ). 

As before, the quantity dh(J^\q,p)) is a legitimate vector held if K pa is sufficiently 
smooth. Moreover, when the Hamiltonian ifC) = h PA o is C 1 we may evolve Hamilton’s 
equations to obtain solutions. 

Proposition 4.5. Let p > | + 2 and a > 0. Then ifC) — h p a 0 is C 1 and given by the 
expression 

H W (q,p) = X ~_P ai K lJ (q a - q b )p bj 

+ pJ% k i{ d kK ij ){q a - q b )p bj 

- \Pa i n Qa l n( d l d kK ij )(Qa ~ Qb^mPb™- 
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Figure 2. (top) Quiver plots of 2D vector fields induced by first order jetlets. 
( bottom ) Streamline plots of 3D vector fields induced by first order jetlets. 


If (q,p)(t) G T*Q$ is a solution to Hamilton’s equations with respect to H^ l \ then J^\(q,p)(t)) G 
£ div (R n )* is a solution of (|4]), and J^\(q,p)(t)) G £ d iv(® n )* is constant in time. 


The proof of the above proposition is identical to that of Proposition 4.2 See again the 


remark following Theorem 4.8 for more details on the kernel smoothness condition on p. 

As before, it is useful to interpret the trajectory q(t) of the previous proposition in 
terms of the curve (ft G SDiff(M n ) obtained by integrating the time-dependent vector field 


dh p ^J^\(q,p)(t )) = K P)lJ * J^\(q,p)(t)). If one chooses po to be the identity or any other 
clement of SDiff(M n ) that satisfies po • (?(0)) = g(0), then q(t) = p t ■ q( 0). This implies in 
particular that the qff\t) are the trajectories of the particles as they are swept along by the 
fluid flow. 

Various vector fields for large kernel smoothness p are depicted in Figure [2] for different 
initial values of the traceless matrix /i/ = pfq\. 


4.4. Ffigher order particles. In this section we will introduce a hierarchy of particle-like 
solutions whose k-th level includes the solutions at the (k — l)-th level. The zeroth level in 
the hierarchy consists of the standard particle-like solutions, while the first level describes 
the particles with internal SL(n) variables discussed in the previous section. The particles 
in the k-th level carry the coefficients of k-th order Taylor expansions of diffeomorphisms, 
or jets. Therefore, we call these particles k-jetlets. 

Let p G SDiff(M n ). The zeroth order Taylor expansion of p about 0 is <£>(0), and the 
collection of such Taylor expansions is all of M n . The first order Taylor expansion of <p about 
0 is 


p l (x) = ^(0) + djp\ 0)x J + o(||x||). 
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The tuple of coefficients (</?(0), Dcp(0)) G M n x SL(n) is called the first order jet of p evaluated 
at 0. Going further, the second order Taylor expansion of p about 0 is 

p l {x) = + djip l ( 0)x j + ^d jk p l (0)x j x k + o(||x|| 2 ). 

We see that D 2 p( 0) is a tensor of rank (1, 2), which is symmetric in the lower indices. We call 
the space of such tensors S%, and we see that the space of second order jets is a submanifold 
Q^ 2> c R n x SL(n) x S^. Finally, for k > 2 the space of k-th order jets is a submanifold 

Q {k) c x SL(n) xSjx^jx-x S l k , 

where Sj is the vector space of (1, j)-tensors which have been symmetrized in the covariant 
indices. The space is equipped with the fiber bundle projection n^ —» M n , which 

projects onto the component. In fact, is a trivial principle bundle where fibers, 
contained within SL(n) x x • • • x Si form a jet group [ DJR13 ]. The jet group serves 
as a finite-dimensional model of the diffeomorphism group (see [ KMS99 , Chapter 4] for a 
description of the group multiplication), and this motivates our interpretation of jetlets as 
models of self-similarity. We define the space for an TV-tuple of fc-jetlets by taking a product 

Q ( N = {(<?!, ■ ■ •, Qn) e Q [k) x • • • x qP I n {k \q a ) * {k \q b ) when a ^ b}. 


(k) 

We coordinatize Q N as follows. We will use Greek indices to represent spatial multi-indices 
on M" (see Appendix C for our multi-index convention). A typical coordinate on will 
therefore look like q a l g where a G {1,...,TV}, i G {l,...,n} and (3 is a multi-index on 
M n . This coordinate is used to model the partial derivative of the z-th coordinate of a 
diffeomorphism at some point z a G M n , i.e. dpp l (z a ). For the statement of the next propo¬ 
sition, recall that the definition of iso(z) was given earlier in ([ 5 ]). Note also that we write 
Jet* : SDiff(M n ) —» Q^\ ] for the function which evaluates the spatial derivatives of a diffeo¬ 
morphism up to order k at the location x G ML. In the subsequent sections we will also use 
the obvious generalization to multiple locations z\,...,Zh G M n , which we will denote by 
Jet* : SDiff(M n ) -A Q {k) . 


Proposition 4.6. The group SDiff(M n ) acts on by a left Lie group action. If z\ ,..., z^ G 
R n are distinct points, then the isotropy group iso(z) C SDiff(M n ) acts on by a right Lie 
group action. 


Proof. Let pi,P 2 G SDiff(M n ). The partial derivatives of (pi o p 2 ) are given by the Faa di 
Bruno formula 


M 


d Q ((fi o p 2 ) = Y^ 


k =1 


\ 


E 


■ ■ ■ d iM k 


\[7lv>7fe]en(a,fc) 




where we wrote II (a,k) for the set of k-th order partitions of a multi-index a. We refer 
to Appendix [C] for the details of our index conventions and to |CS96, Jacl4j for a precise 
description of the multivariate Faa di Bruno formula. One can read off from the expression 
that a fc-th order derivative only depends on /c-th and lower order partial derivatives of p\ 
and p> 2 - A left SDiff(M n ) action is induced on by setting p ■ q = Jet*(<^ o fj) for any if 
such that Jet*(i/T) = q. That this is independent of the choice of if follows from observing 
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that the Faa cli Bruno formula only uses data in q = Jet k (ip) and nothing more. In local 
coordinates, this action takes the form 


(<p ■ q)' l a = 


M 

E 

k =1 


/ 


\[7i>-»7fc]en(a,fe) 


d, 


31-"3k 




(q ,0, k 


,31 


71 


Q Jk 
d 7 k 


except for the component where |a| =0 in which case we observe (p ■ q) 1 = p l (q^). By the 
same construction, a right iso(z) action is induced on by setting q ■ p — Jet k z {iJ) o p) for 
any i/j such that Jet^(V’) = <?■ We can choose distinct points ..., G R n and apply the 
same process to ■ In this case we observe the action to be 


(q ■ ¥>)“ 


M 

E 


/ 


v r Uu . 

\[7i,...,7fc]en(a,fe) 


\ 

,j k ] d 'n ( f jl ( Z a)--- d ' r y k (Za) 


□ 


Just as in the previous sections, the actions of SDiff(M n ) and iso(z), which commute, lift to 
actions on T*Q^\ The associated momentum maps are given implicitly by how they act on 
the respective Lie algebras. In particular, the left action of SDiff(R n ) yields the momentum 
map 


( J ( L k) (q,p),u) = 


N 

E 

a= 1 


/ 


PamU m (qj? } )+ Y P o,m 

\a\<k 

\ 1<£<|o| [7i,...,7^]en(a/) 


'd. 


31 — 3e 


u 




31 


71 


q • le 

Ha r/t 


for arbitrary divergence free vector fields u G 3Cdiv(®i n )- Equivalently, we may define J 4 as 
the unique map such that 

(10) (J^\q,p),u) = (( q,p),Jet k z (uop )) 

for any p G SDiff(M n ) whose /e-jet is given by q for any (q,p) G T*Q$. 

The right action of iso(z) yields the momentum map defined implicitly by the relation 

N l / 

(4 fc) (?,P), w) = Y Y Y Pai a qa[n-M Y a 7 m ™ jm ( Z a) J] ^ 

a=l |a|<fc m=l \n^m 



where w G 3Cdiv(®i n ) is such that w(z) = 0 (this describes the Lie algebra of iso(z)), and 
is the natural generalization of the Kronecker-delta symbol to multi-indices. Equivalently, 
we may define as the unique map such that 

(11) (Jft\q,p),w) = ({q,p),3et k (Tp-w)), 


for any p G SDiff(M") whose /e-jet is given by q for any (q,p) G T*Q$. Here, we wrote 
Tp ■ w for the function obtained by applying the differential of p to w. 
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(k) (k) 

We can write J L and J R explicitly, using the Dirac-delta distribution, as 

N ( 

Jl\q,P .) = ^2 I Pamdx m ® 5 q (0) 

a =1 \ 

+ V q J\ 

|o|<fc 

l<i<\a\ [ 7 i,..., 7 £]en (a,£) 

and 

N l / 

(12) 4\ q ,p )-V V 5] pav*]E'- 1 )' 7 " 1 n J:' 

a=l |a|</c m=l \n^m 

Proposition 4.7. The maps J2 and J^ form a weak dual pair. 

Proof. The result is a direct application of (GBV12I Corollary 2.8]. To make the exposition 
more self-contained, we provide some details. By ] CB\ 1 2. Corollary 2.6] we need only 
show that Jl and Jr are equivariant, and that Jl is invariant under the right action of 
iso (z). Equivariance follows from the fact that Jr and Jl are derived from cotangent lifted 
group actions | IAM08. Corollary 4.2.11], So we need only illustrate that Jl is (right) iso(^) 
invariant. This can be seen as a consequence of the commutativity of the left and right actions 
on Q^ k \ For notational clarity, let us denote this right action by p : Q ® x iso(z) —» Q^ k \ 
Explicitly, any element q G is expressible as the fc-jet of some diffeomorphism p, and 
p(g,-0) = Jet^cp ° ip). 

Similarly, any element of the tangent fiber T. f SDiff(M n ) may be written as a composition 
u o p for some u G 3£div(BC)- Elements of TQ^ are of the form Jet^(w o tp) for u G Xdi v (®C) 
and p G SDiff(M n ). Given this representation, the tangent lift of the action p, also denoted 
p : TQW x iso (z) —> TQ^ k \ is given by 

(13) p((q, v), ip) = Jet k (u op oif) 

for (q,v) G TQ ^ and where u G Adiv(BC) and p G SDiff(M n ) are arbitrary up to the 
constraint (q,v) = Jet k [u o p). The cotangent lifted action is a left action, p* : iso(z) x 
T*qJ) —> T*Q^ k \ defined implicitly by the condition 

( 14 ) (p*{ip, (g,p)),p((g,u)^ -1 )) = ({q,p),(q,v)) 

for all (q,p) G T*Q^ and (q,v) G TQ^ k \ This action is equivalent to the one defined in the 
discussion preceding [ AM081 Corollary 4.2.11], In particular, p*(ip, (g,p)) is a covector over 
the point q = p(q, ip^ 1 ). By the definition of Jl we observe 


d ^ 7 




... A M f) . 
i lci ~u U n- 


■3r 


Qa 


( 0 ) 


(J L {p*(ip, (q,p))),u) = (p*(ip,(q,p)),Jet k (uop)), 
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where p G SDiff(M n ) is any diffeomorphism such that Jet^(<^) = p(g, "0 1 )- If we let ft be 
such that Jet^(<£>) = q then we can simply choose ft = cp o ft~ l . Thus we get 

(J L (p*(*P,(q,p))),u) = (q,p)), Jet^uopoft- 1 )) 

= (p*Wh (9,P)),p(Jet^(uo^),'0“ 1 )) 

= {(q,p), Jet*(uoy>)) 

= (JL(q,p),u), 


using (13) and (14) in the second and third equality. Thus, we see that Jl is invariant under 
the right action of iso(z) on Q^ k > and the result follows. □ 


In the case of k — 1 we obtain the weak dual pair of the previous section and Proposition 
4.4| is a corollary of Proposition 4.7| Proposition T7 gives us the final result on jetlet 


parametrized solutions. 


Theorem 4.8. Let p > | + k + 1 and a > 0. Then = h. [K(J o ./[^ is C ] . Let x(t) 
be a solution to Hamilton's equations on T*Q^\ then J^\x(t)) is a solution to Hamilton's 
equations on 3C div (M n )* an d is constant in time. 


Remark that K pa has smoothness C 2p ~ n ~ l by considering that the Fourier representations 
of its derivatives are integrable, while we need K pa G C 2k+1 to allow composition with k-th 
derivatives of delta distributions on each side and still obtain a C 1 Hamiltonian. 

As before, an even richer family of velocity fields is generated by a single particle of this 
type. At order k — 2 this yields four new varieties of velocity fields per particle. Two 
examples of such velocity fields are depicted in Figure [3| 



-4 -2 0 2 4 



Figure 3. Some velocity fields generated by particles of order k = 2. 

4.5. Kelvin’s circulation theorem. In this section we relate the conserved quantities 
associated with to Kelvin’s circulation theorem. Let ip t G SDiff(M n ) denote the flow- 
map produced by a solution to Euler’s equation. Let 7 (s,t) = 77 ( 7 o( s )) f° r some loop 7 o(s). 
Kelvin’s circulation theorem states that the circulation 

r (t) = j> u(^(s,t)) ■ d s 'f(s,t)ds 

is constant in time. It was shown in [Arn66| that this conservation law is an instance of 
Noether’s theorem. In particular, circulation is one of the conserved momenta associated with 
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the particle relabeling symmetry of fluids. More specifically, recall that ho is the Hamiltonian 
for Euler’s equations, and this Hamiltonian is SDiff(M n ) invariant. Moreover, we have a 
weaf0 dual pair J spat i a i : T*SDiff(M n ) —> Xdi v (^ n )* and J conv : T*SDiff(M n ) —>- iso(z)* of 
spatial and convective momentum maps induced by the left action of SDiff(M n ) and the right 
action of iso(^) C SDiff(M n ) on SDiff(M n ) itself. By applying Theorem A.9 with A = J spatial 
and J 2 = Jconv we know that Hamiltonian dynamics on T*SDiff(M n ) with respect to a 
Hamiltonian of the form ho J 1 will exhibit the constant of motion 

(15) (^conv(^,P v ),w) = {Pv,T(p ■ w) for all w G \so(z). 

This instance of Noether’s theorem applies to regularized models as well as to the non- 
regularized case of Euler’s fluid equations. In the case of Euler’s fluid equations, the equiva¬ 
lence between the above conservation law and Kelvin’s circulation theorem is demonstrated 
by a heuristic argument jAK98 . Chapter 1, Theorem 5.5], which goes as follows. Consider a 
given closed curve 70 : S 1 —* M n and a family of vector fields w £ such that w £ ( n f 0 (s)) = 7 p (s) 
and such that the (weak in L 2 ) limit w = lim e _ ) .o' u; £ is the generalized function 

w(x ) = 7 q (s) <g> 5(x - 70 (s)). 

Assuming £ is small and writing p v = (<p, ■ )x, 2 , we hnd 
{Jconvi^^Pi P(p) -i 'We') {Pi Tp ‘ W £ )l,2 

= / {ip(x),Dip(x) ■w £ (x)) VL ndx 


(0(7o(s)), Dy>(7o(s)) • 7 o( s ))k- ds 


= cb{u('y(s,t)),d 8 'y(s,t)) Rn ds, 


where u = ip o Lp~ x is the usual Euler representation of the fluid flow and 7 (s, t) := p t (p'o(s)). 
Therefore, conservation of J con v leads to conservation of circulation. 

A more rigorous correspondence is developed in [HMR98] . Any m G Xdi V (K n )* can be 
written as a one-form density fh ® p where /i is the volume form on M n . For any smooth 
curve 7 : S 1 —> M n we may consider the current J^(y) G Xdi v (K n )** defined by 

(J ^(7 ),m) = rh. 


By Theorem 6.2 of |HMR98j , if 'yt(s) = <ft( 7 o(s)) and m = fh ® p satishes the ideal incom¬ 
pressible fluid equation (perhaps regularized), then the circulation 

T(f) = (J^(y t ),m(t)) 

is constant in time. It is in this sense that Kelvin’s circulation theorem follows from the 
particle relabeling symmetry for any SDiff(M n ) invariant Lagrangian. 


In light of this discussion, it is natural to consider (15) as the fundamental conservation 


law. The main goal in the remainder of this section is to show that the jet-particle solutions 
satisfy conservation laws that are ‘shadows’ of this fundamental law in the sense that they 
are associated with a partial relabeling symmetry. To see this, it is useful to define the 
k -th order isotropy group, iso ^(z) = {-0 G SDiff(M n ) | Jet^(0) = Jet^(id)}, where we 


5 With SDiff(R”) as right action this would be a proper dual pair, but isofz) C SDiff(R ra ) does not act 
transitively on the level sets of J spa tiaii hence the dual pair is weak. 
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wrote id for the identity mapping and z is shorthand for z±,..., Zn £ M n . Then we see that 
= SDiff(M n )/iso^(z), and the corresponding quotient map is given by Jet*. Let us also 
introduce the right action R : SDiff (M n ) x iso(z) —> SDiff (M n ) given by right composition of 
functions and write so that R^p = p o - 0 . We also introduce the cotangent 

lifted right action TR* by means of the defining relation 

(16) (T 'Ry-ip v , v v o ij>) = fa, V{p ) 

for arbitrary p v £ T* SDiff (R n ) and v v £ T v SDiff (M n ). 

The space T*Q$ naturally embeds into T* SDiff (M n )/iso^(z), where the quotient is by 
the cotangent lifted action. More precisely, for any ( q,p ) £ T*Q$ we can construct the 
corresponding element i(q,p ) £ T* SDiff (M n )/ iso^(z) in the following way: take any p such 
that Jet*(</?) = q and find L,p(q,p) £ T* SDiff(R n ) that satishes, for all v v £ T v SDiff(M n ), 

( 17 ) Mq,p), v v ) = ({q,p),{v v op- 1 ) -q), 

cf. [MMO + O7l Equation (2.2.4)], where we denote by (ry, o p _1 ) ■ q the infinitesimal action 
from the left of o p ~ 1 £ Xdi v (K n ) on q. Then set i(q,p) = [< v (g,p)]. 

To see that i(q,p) is well defined, note that if p 1 = p o ib for tb £ iso ^ k Hz), then for any 
V e T P SDiff (U 1 ) we have 


(TR%,-iL v (q,p), vp) = (‘v(?)P)iV 0 ^ X ) 

= ((q,P), (vp o ij;- 1 o p- 1 ) ■ q) = {(q,p), (vp ° p'~ l ) ■ q) = (tp(q,p),Vp) 

using (16) and 0 in the first and second equalities. Since [TR^ 1 t lf (q,p)\ = [tp(q,p)\, we 
conclude that i(q,p ) is well defined. 

We claim that 


(18) i(q,p) = [TRpp\q,p)}, 

where p is such that Jet k z (p) = q. This follows since for any v r £ T v SDiff (R n ) 

(Ti?;_ijf(g,p),^) = (JL ) (q,p),v v o(p- 1 ) = ((q,p),(v v op~ 1 ) ■ q). 

With these preliminary remarks in mind, we can now show the commutativity of the 
following diagram: 



Here, S — (n 1 o i) (T*Q^^ and J spatial and J con v are defined in the natural manner through 
the left and right actions of SDiff (W 1 ) and iso(z). To verify the left side of the diagram, we 
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let (q,p) G T*Q$ and p such that Jet*(<£>) 
weX div (R n ) 


q. Then we obtain from (18) that for any 


([Jspatial](i(g,p)),w) = {J B p a tAal(TI? v -ij¥ ) (q,p)),w) 

= (TR* v - 1 J i L k) (q,p),wo(p) = (jP(q,p),w), 

as required. This also shows that [J spa tiai] is well defined on T* SDiff(R n )/ iso^ k \z) since the 
final expression does not depend on the choice of representative p for q. 

For the right leg, note that 

<[^conv](i(?,P)),w) = <Jconv(Ti?;_iji%,p)),U;> 

= (TR*^ J {k \q,p), T(p ■ w) — (J { L k \q jP ), ( Tip • w) o <p -1 > 

= {(q,p),Jet k z (Tp-w)) = (Jft\q,p), w). 


Here we used (18), the definition of J con v, and ( |16| ), respectively, in the first three equalities, 
and ([TTJ) in the final step. Again we see that [Jconv] is well-defined since it does not depend 
on p. More explicitly, for another representative p' — p o 'tjj with 0 G iso^(z) we see that 
-0 gets projected out by (the tangent map of) 7r = Jet^. That is, the restricted momentum 
map J conv | s is invariant under the right action of iso^(z) and therefore descends to a map 
[Jconvls] defined on i(T*Q\$). 

The right legs of both weak dual pairs yield the conserved quantities. From the diagram 


it is clear that the conservation of J 


(k) 


R 


conservation of J conv in (15). Since J c , 


exhibited in our particle models is a shadow of the 
nv corresponds by Noether’s theorem to the (large) 


subgroup iso(z) of the right symmetry SDiff(M n ) that generates conservation of circulation, 

(k) 

we see that conservation of J n ; in the jetlet solutions is a shadow of the conservation of 
circulation. In other words, our particle models contain a model of Kelvin’s circulation 
theorem within them (cf. [ D.TR13 , Theorem 5.5]). 

This diagram also provides some insight into the relationship between the developments 
in this paper and classical Marsden-Weinstein reduction theory [ MW74] . For instance, 
letting Jconv : T*SDiff(M n ) —» iso^(z)* be the momentum map associated to the cotan¬ 
gent lift of the right action of iso ^(z) on SDiff(R n ), one can show that i(T*Q = 

(Jconv ) -1 (0)/ iso^(2). For more details on the connections between general reduction theory 
and the results of the present paper, see Appendix [Bj 


5. Particle mergers 

In this section we discuss some explicit dynamical behavior of the particle model, in 
particular we study ‘collisions’ of jetlcts. For the zeroth order particles case, this was already 
analyzed by Mumford and Michor |MM13bj . and they found that two particles can merge 
in infinite time, or bounce off each other, depending on the ratio of their relative angular 
and linear momenta. To find the explicit behavior analytically, we shall restrict to two 
dimensional space and two 0-jetlet particles with zero total linear momentum. We identify 
the asymptotics of the merged state as the dynamics of a single 1-jetlet particle. 
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We start with the Hamiltonian H *+ for two 0-jetlet particles, 


H = A Pai Kt3 (Qa - %)Pi 


bj 


a,b= 1 


on T*M 2n with the canonical Poisson brackets. Translation symmetry allows us perform 
symplectic reduction. We switch to a center of mass frame by choosing new coordinates 


(19) 


q = 2(91 + 92), 


q = q2-qi 


with canonically associated momenta 

P = Pi + P 2 , 


P = 2(^2 - Pi)- 


The Hamiltonian in these coordinates becomes 

H = ^(A^'(O) + (q))pj + pi( (0) - K^(q))p r 

When n — 2 and the total momentum is zero, i.e. p = 0, we can perform another symplectic 
reduction by rotational symmetry. We switch to polar coordinates 

(20) q = (r cos(0), r sin(0)) 

with canonically associated momenta 


pT _ (' COS(0) 


sin ( 0 ) 


Pr 

Pt 


— Rd 


,sin(^>) yM 

where is a rotation matrix. In these coordinates the Hamiltonian is given by 
H = p [K ! ’fm - K l ‘(q)}p, = (ilb Rl[K( 0) - A'ifl, . (r,0))]fl* (gj 

= (|j [A(0) - K(r, 0)] (gY 

where in the last step we used that K as a tensor is invariant under rotations. Since 0 is a 
cyclic variable, we find that its associated momentum 


p^ = —r sm((j))pi + r cos (0)p 2 = q Ap 

is conserved. Remark that this relative angular momentum is not exactly the total angular 
momentum when p ^ 0. 

Let us now choose the smooth kernel K = iL 0O; i given (up to a scaling factor) by 

K"\x) = (e-- - T(i - e -P)y + (1(1 - e-o) - e“-) 

II II2 

where p := Note that AT y '(0) = p,S ij and {{)) = 0. Using rotational symmetry we 
set 0 = 0, and obtain 

(21) B = f(l - -“■')) + So-2 e ‘' + ^(l-«■'))• 

This system is Hamiltonian in (r,p r ) with p$ a parameter, so the level sets of H determine 
the motion. For p<f> = 0 we see that H oc p 2 r r 2 /(r 2 ) where / is an analytic function with 
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/(0) > 0. Thus the level sets of H near r = 0 look like hyperbola with p r = 0 and r = 0 
as asymptotic axes, see the left image in Figure Q. Hence, two particles approaching each 
other head-on will ‘collide’ in infinite time; even though their momentum blows up, their 
relative velocity decays exponentially. 

If P<t> 7^ 0 then the first non-constant contributing term has sign opposite to the terms 
involving p r , hence we see a new region being created in the right image of Figure (JdJ) where 
orbits approach r = 0 and then ‘bounce off’. There are two asymptotic orbits that approach 
the r = 0 axis at finite momentum; the limit point can be calculated from the fact that 
p r /r = 0 must hold there in the limit r —» 0. We find that p r = ±^/5j6p<p. From Figure [d] 
it is clear that this is a good approximation for the asymptotic curve. 




FIGURE 4. Contour plots of the reduced Hamiltonian (21) with r horizontal 
and p r vertical, left for = 0 and right for p= 0.3. The red line shows the 
positive asymptotic value for merging or bouncing. 


Finally, we can reconstruct the asymptotic 1-jet let trajectory that these two merging 0- 
jetlets converge to and verify that this trajectory is indeed a solution of the Hamiltonian 
vector field for a 1-jet let. 

To analyze the asymptotic behavior, we expand H around p = 0: 

Since H and p 2 ^ are preserved, we can solve for p r in terms of r, and we find 

2 PPl = (8 H - 3pj) (1 + 0{p)), 
rp r = \J&H — 3p^ + 0(r). 

We write ( := yJsH — 3 p 2 ^ and thus obtain asymptotically p r = £ + 0(1). Further, we have 
dynamics 

dH 1 2 4 . • dH 3 

r= 9U = 4 Prr ' + ° (r) ' 4 ’ = W t = l P ^° i) 
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and reconstruct 

9 = 0, ^ =rR *'\n)’ 

p = 0, = =^'(pj+°(i)- 

Now we consider the image under Jl of the asymptotic solution curve. By testing against a 
vector field in 3C div (M n ) we find 

JL(qi,Pi,q2,P2) =pi ®S qi +P2®^2 

= p®8q+p® (-tfdjSq) +o({\p\ + |p|)|g| 2 ) 

and inserting our reconstructed solution leads to 


( 22 ) 




i 


dx 1 dj5g + 0(r) 


with (j>(t ) = | p,/, t. Our aim is to show that the factor in front of dx 1 djS, which we shall 
denote by p/(i), corresponds to p/(f) for a 1-jet let with position q ^ 0) = 0 and momentum 
= 0. From here on, we use the Frobenius inner product to identify 51(2)* = 51(2), which 
for explicit matrices corresponds to taking the transpose. For such a setup we have equations 
of motion 


Pi = Pi k d k u*( 0) - <V(0)/V = -d km IV\d)^rpi + d im K kl ( 0)/W', 

or in short, p = — ad^(p) = [/x, £ T ] with = djU l { 0) = —djhK tl ( 0)/x/\ To ease calculations, 

let us choose the basis 


CO = 





and note that these matrices have norm \/2. We decompose /i = j co + s cr+tr and the second 
derivative of the kernel as a tensor product. A calculation verified by symbolic computer 
algebra software shows that 


dijK kl ( 0) 


/ 


1 

4 


V 



o 1 
1 o 

3 0 
0 1 


\ 


1 k i 


^ / k l j k—l\ 

A \ a i a j T i T j ) ’ 


where the indices i , j and k,l label the outer and inner matrix elements respectively. With 
these decompositions we find that £ as a function of /; can be written as 


£ = j u + 


1 

2 


(sa + tr). 


Using the commutation relations [c o, a] = 2r, [a;, r] = —2cr, [a, r] = —2a;, it then follows that 

1 3 

(23) p = [/x,£ T ] = [ju + sa + tr, -j u + -(sa + tr )] = -3 j(-ta + ST) = --j [aa,/x]. 


That is, the (a, r) components of p as a tensor rotate with angular velocity —3 j. 
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On the other hand, from (]22|) we have for the asymptotic dynamics that 

c 


Kt) = 


Ra 


p<t> 


Ra 


= - Ad 


Ra, 


c 

Pd 


with only 0 = |depending on time. Note that (1 st(2) as a matrix, but since it is 

actually a dual element, we can simply ignore its trace part and project it out. Also note 
that the u> component of J2(0) is j = \p Differentiating with respect to time yields 

. . _ 3 3 

P = -0 ad u (p) = /i], 


and comparing to (23) we find that the asymptotic solution of the two merging 0-jetlets 


matches that of a 1-jetlet with the same angular (and linear) momentum. 

Let us finally suggest a more abstract way to view these particle mergers. Our hierarchy 
of reduced spaces T*Q$ embeds into 3Cdiv(®0)* under the momentum map Jl- Consider 
a merging pair of 0-jetlets, described by a curve x 0 (t) G T*Q^ ■ As the particles approach 
each other, x 0 (t) approaches the boundary of T*Q^ given by 

d{T*Q < 0) ) = {(qi,Pi,q 2 ,P 2 ) G T*R n x T*R n \ q x = q 2 }. 

On the other hand, the image curve yo(t) = Jl(^o(^)) G ^div(^ n )* consists of two covector¬ 
valued delta distributions at qi,q 2 , and in the limit as their distance goes to zero, this can be 


approximated by a momentum valued distribution of a delta and its derivative (22), that is, 


an element yi(t) = Ji(xi(f)), where X\(t) G T*Q^ is a curve in the space of single 1-jetlet 
particles. 


* div K n r 


2 °) as a subset of T*Q± , and consider a topology on 
in which the embedding is continuous, see diagram (24). This picture naturally 


We can view the boundary of T*Q 2 


(k) 

generalizes to the whole hierarchy of spaces T*Q N , suggesting that it might be interpreted 
as a CW-complex. In this setting, the question whether the solution curve yo(t) of the 
merging 0-jetlets converges to a solution curve yi(t) of a 1-jetlet, basicall}^] boils down to the 
question whether the vector held of the dynamics on Xdi v (®D)* is continuous. We have not 
pursued in detail the question of which topology on 3Cdi v (^ n )* t° use for this more abstract 
characterization. 


(24) 



6 One has to be careful, however, since continuity of the vector field will only imply that the asymptotic 
curve yo(t) £ Jl(T*Q^) is a pseudo orbit of the 1-jetlet dynamics. This does not imply existence of 
a solution curve yi(t) £ Jl{T*Q j 1 -*) that yo{t) is asymptotic to; that would require the ‘limit shadowing 
property’, which is closely related to hyperbolic properties of the dynamics [PPT121 IRibl4j . 
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6. Numerical experiments 

We have performed a number of numerical simulation^] of jetlet particles. These confirm 
that the conserved quantities are indeed preserved and that two 0-jetlet particles merge as 
shown by the analysis in Section [5] providing a sanity check for the formulas in the previous 
section and the numerical code. Moreover, the merging behavior shows to be stable under 
perturbations of initial conditions. 

The numerical code works for any spatial dimension n > 2, but for the sake of tractability 
and simplicity we have studied n = 2 

As basic experiment we take two 0-jetlet particles with initial states 

Qi = (-3,0), pi = (1.5, —d), 

92 = (3,0), P 2 = (—1-5, d), 


aimed at each other with an offset parameter d. The initial state is given in center of mass 
polar coordinates, see (19) and (20), by 


r = 6, p r = —1.5, 

0 = 0, p<f, = 6 d. 


Furthermore, we use a — 1 throughout our experiments. The experiments show that for 
d = 0.27 the two particles merge while spinning around each other, while for d = 0.288 they 
get close, but then emerge from their close spinning state and scatter in opposite directions. 
This confirms the analytical value of d = ^/3/40 ~ 0.2739 within reasonable precision, 
noting that this is the asymptotic value for particles starting close to each other. 

We performed a number of more complex simulations, all of those being small perturba¬ 
tions of the basic experiment described above. First, we added a small angular momentum 
‘spin’ component to both particles, turning both into 1-jetlets, one level higher in the hier¬ 
archy. Then we added a third particle (both a 0-jetlet and 1-jetlet) at such a distance and 
momentum that it exhibits medium range interaction with the first two particles. Finally, 
we added a small hyperbolic-like ‘stretching’ momentum to the first particle only. We found 
an analytic study of these configurations to be infeasible, but the simulations show that the 
behavior observed in the basic experiment persists. We can find parameter values of d close 
to the original one where the system shows a transition between the two particles merging 
or scattering. 

These experiments also confirm the preservation of the conserved quantities present in 
the system. For all experiments described above we observed that the energy, total linear 
and angular momentum, as well as Jr (individually for each particle) were preserved with 
absolute errors less than 4 • 10~ 4 over a time of 60 seconds, while the energy was of the 
order one. Unlike Jr, which is conserved for each particle, linear and angular momentum 
can be exchanged between particles (although total momentum is conserved). Figure [5] left 
shows the angular momentum of two scattering 1-jetlets and the right plot of two scattering 
0-jetlets interacting weakly with a third jetlet particle. 


7 The simulation was written using Python and NumPy, the source code and generated videos can be 
found at: https://github.com/hoj 201/incompressible_j et_particles 
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Figure 5. Angular momentum exchange between two scattering 1-jetlets 
(left) and with a third particle present (right). 

7. Conclusion 

In this paper we derived a hierarchy of weak dual pairs which induce a family of particle-like 
solutions, called jetlets, and conserved quantities that shadow the conservation of circulation. 
The jetlets have internal degrees of freedom given by a jet group. As the jet group is a finite¬ 
dimensional model of the diffeomorphism group, we suggested the use of jetlets as a finite¬ 
dimensional model of self-similarity, wherein a “large” diffeomorphism advects a “small” 
diffeomorphism. We also studied the dynamics of mergers and provided a rigorous analysis 
showing that merging 0-jetlets asymptotically approach 1-jetlets. 

The developments discussed in the present paper give rise to a number of promising 
directions for future research. These include: 

(1) An investigation of the relationship between jetlets and point vortices or vortex blobs. 

(2) Further investigation of the numerical implementation. The use of parallelization 
and the fast multipole method would be particularly interesting to consider. 

(3) Finding a way to implement boundary conditions. In such scenarios, the kernel is no 
longer invariant under rigid transformations and we must consider a general kernel 
K : M x M —> M nXn where M C 1" is an rt-manifold with boundary. 

(4) An analysis of convergence to Euler equations when a —> 0 for the case where the 
power of the Helmholtz operator, p, goes to infinity. The advantage of the p = oo case 
is that the limiting kernel can be written in terms of elementary functions [ MGl4j . 
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Appendix A. Hamiltonian mechanics 


The goal of this section is to prove Theorem A.9 (see page 31). Those who understand 


and accept these theorems on a first reading should be able to skip this section without any 
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consequence. Most of this section will be a crash course in Poisson geometry and Hamiltonian 
mechanics as described in [ AM0811MR991 IWei83j . 

A typical introduction to Poisson structures in mechanics begins by considering Hamilton’s 
equations 

. dH . dH 

q ~~&p ’ 

If we consider the two-form dpAdq, then Hamilton’s equations can be written as (q,p) _I (dp A 
dq) = dH(q,p). This is the starting point for symplectic geometry, which will be discussed in 
Section A.l Alternatively, these equations can be written using the bilinear map {•, -} can : 
C\T*Q ) x C 1 (T*Q) -A C°(T*Q) given by 

_dF8G _dG dF 
tr j- can — — — — . 

dq op dq op 

In particular, we may write q = {q, H} can and p = {p,H} can . The object {•,-} can is a 
special case of more general object known as a Poisson bracket which will be introduced in 
Section IA.21 

A word of warning: symplectic geometry has developed greatly since its origins in me¬ 
chanics, and has branched into an independent subfield of pure mathematics. Many notions 
were revised and optimized in the 1970’s and 1980’s for the purpose of proving theorems. 
Occasionally these revisions entailed a sacrifice in clarity, from the perspective of “outsiders”. 
This paper is intended to allow “outsiders” (such as ourselves) to reap the benefits of Poisson 
geometry. Therefore, we will cut away as much abstraction as possible in this introductory 
section. Nonetheless, a minimal amount of abstraction is needed in order to maintain math¬ 
ematical rigor and stand firmly upon the shoulders of giants. 


A.l. Symplectic manifolds. We begin with the definition. 

Definition A.l. Let S be a manifold and let uj be a closed two-form on S such that the map 
“v £ TS i-A u(v,-) £ T*S” is weakly non- degenerate^ We call u ; a symplectic form. We 
call the pair ( S,u >) a symplectic manifold. 

All of the expressions derived in this article are formal, and we refer to GBV 121 for the 
functional analytic details of infinite-dimensional symplectic manifolds. As a first example, 
consider the manifold R 2 with coordinates ( q,p ). The two-form dqAdp is a symplectic form. 
Given a manifold Q, the cotangent bundle T*Q has local fiber bundle coordinates given by 
(q 1 ,..., q n ,Pi ,... ,p n ) and there is a unique symplectic form which is locally expressed by 
dp 1 A dq . t , where a sum on repeated indices is assumed. This local expression corresponds 
to a global symplectic form on T*Q , known as the canonical symplectic form and denoted 
cu ca n [ AM08I . Theorem 3.2.10]. In fact, given any symplectic manifold (S,oj), the dimension 

of S is even, and there exist local coordinates (q 1 ,..., q n ,Pi, ■ ■ ■ ,p n ) such that ui lo = ly dpi A 
dq 1 . This is known as Darboux’s theorem and we call this type of coordinates Darboux 
coordinates |AM08l Theorem 3.2.2], 

Given a function H : S —> R, the exterior derivative is the one-form dH : S -A T*S 
expressed in local coordinates by dH(x ) = |^dx*. The Hamiltonian vector field X H : S -A 
TS is the unique vector field defined by the condition X H _lu = dH. The symbol “ _J ” is 
the operation of contraction between the contravariant indices of Xh and the first set of 

8 A linear map L : V —> V* is weakly non-degenerate if L is injective. If V is finite-dimensional, this 
simply means that L is invertible. 
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covariant indices of u. In Darboux coordinates, the Hamiltonian vector field induces the 
equations of motion q l = = — |p. 

An important aspect of study in Hamiltonian mechanics is that of symmetry. This yields 
the following notions. 

Definition A.2. Let G be a Lie group and let p : G —>• Diff(S') be a group action on a 
symplectic manifold (S', cu). The group G is said to act symplectically if 

v(p(g) -v,p(g) -w) =u(v,w) 

for any g G G,v,w G T X S and x G S. If g is the Lie algebra of such a group, the momentum 
map, J : S —> g*, is defined by the property 

d(J,£) =£,s- iw. 

Alternatively, we can characterize a momentum map J : S -A g*, as the unique map such 
that X(j£) = fs for any (eg. In the special case where S = T*Q , a left/right action of G 
on Q can be lifted to a right/lcft symplectic action on T*Q given by 

(q,p)eT*Q^(g- 1 - q ,g*p) eT *Q. 

where g*p is the unique covector such that (g*p, v) = (p,Tg ■ v). In this case the momentum 
map is characterized by the condition 

(26) (J{q,p),0 = (p,f • q)- 

This is contained in Theorem 12.1.4 of [ MR99] . 

Finally, given two functions /, h G C°°(S) we can consider the function {/, h} = oa(Xf, X h ). 
In Darboux coordinates {/, h] = §^i§jf — Hamilton’s equations can then be written 

as q l = {q\ h}, p t = {p,, h}. We call {•, •} a Poisson bracket, and it is the subject of the next 
subsection. 

A.2. Poisson manifolds. We begin with the definition. 

Definition A.3. Let P be a manifold , and {•, •} be a bilinear operation on C°°(P) such that 
(C°°(P), {•,•}) is a Lie algebra and {-,h} has the derivation property for any h G C°°(P). 
That is to say 

{gf, h} = {/, h} ■ g + {g, h} ■ /, 

for any f,g,h G C°°(P). We call •} a Poisson bracket, and we call the pair (P, {•, •}) a 
Poisson manifold. 

The most important example of a Poisson bracket is that of a Poisson bracket on a sym¬ 
plectic manifold (S,oj). Here the Poisson bracket is {/, g} = uj(Xf,X g ). When S' is a 
cotangent bundle, and u is the canonical symplectic form, we call this bracket the canonical 
Poisson bracket. 

The second most important example of a Poisson bracket, after the canonical Poisson 
bracket, is the Lie-Poisson bracket. Let g be a Lie algebra and let g* denote its dual. The 
Lie-Poisson bracket on g* is given by 

(27) {f,g} s *(x) = ±/x, 

where (•, •) is the canonical pairing between dual-vectors and vectors, and [•,•] is the Lie 
bracket on g. The “+” Poisson bracket is nothing but the canonical Poisson bracket on T*G, 
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mapped to the space g* via the left trivialization map A : (g,p) G T*G i—» ( L g )*p G g*. The 
” bracket is obtained through the right trivialization map p : ( g,p ) G T*G i—> ( R g )*p G g*. 

On a Poisson manifold (P, {•,•}) the derivation property implies that the functional oper¬ 
ator {-,h} is equivalent to the Lie derivative operator of a unique vector held X h : P -A TP. 
That is to say, X h is the unique vector held such that Cx k [./’] = {/, h} for any / G C 1 (P). 
We call Xh the Hamiltonian vector field, and the ODE x = Xh(x) is called a Hamiltonian 
equation. It is standard to write this ODE as “x = {x,H}”, despite the fact that one typi¬ 
cally intends for “x” to represent a point in P, and not a function. Since one can take “x” 
to be a place-holder for a set of local coordinate functions which determine x uniquely, this 
sloppiness is usually harmless. 


Proposition A.4 (Proposition 10.2.2 |MR99j ). Let (P, {-,-}) be a Poisson manifold. Then 

Corollary A.5. Let (S', a;) be a symplectic manifold and let h, f G C°°(S). Then [A fi,Xf\ = 
~ X u(x: h ,x f )- 

Proof. X u ( Xh ,x f ) = x {hj} = ~[X h ,X f \. □ 

Definition A. 6 . Let (Pi, {•, -}i) and (P 2 , {-, -} 2 ) be Poisson 7nanifolds. A map if : P\ —> P 2 
is called a Poisson map if{f°if,g° = {/, < 7)2 0 V’ f or an V L 9 £ C°°(P 2 ). 

Proposition A.7 (Lemma 1.2 of |Wei83j or Proposition 10.3.2 of [MR.99J ). Let if : P\ —> P 2 
be a Poisson map. Let h 2 G C 1 (P 2 ). If x(t) G Pi is a solution to Hamilton’s equations 
with respect to hi — h 2 o if G C 1 (Pi), then y(t ) = fi>(x{f)) G P 2 is a solution to Hamilton’s 
equations with respect to h 2 . 


Remark that unlike in (MR99 ], we only require C 1 smoothness since we do not use existence 
and uniqueness of solutions. 

When the dimension of P 2 is larger than that of P 1; Proposition A.7 allows one to find 
solutions of Hamiltonian equations on P 2 by solving lower-dimensional Hamiltonian equations 
on Pi. 


A.3. Weak dual pairs. In this section we review the notion of weak dual pairs [ GBV12 j. 
This is a relaxation of the more frequently invoked notion of a dual pair [ MW83 . !Wei83] . 
Let (S', a;) be a symplectic manifold. Given a distribution V C TS, denote the hber over 
x G S by V x C T X S. The syinplectic orthogonal to V is the distribution 

V u = {we TS | w(w,v) = 0, Vn G V}. 

Definition A.8 (Weak dual pair |GBV12j ). Let J \: S —> Pi be a Poisson map. The kernel 
of Ji is the distribution 

kernel(Ji) = {v G TS \ TJ X ■ v = 0}. 

If J 2 : S —>■ P 2 is a Poisson map as well, and 

kernel (J 2 ) w C kernel (A) , kernel {Ji) w C kernel (J 2 ) 

we call the diagram 

Pi A S A P 2 

a weak dual pair. 


We would have a proper dual pair if the kernel inclusions were replaced by equalities. 
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Theorem A.9. Let Ji, J 2 : S —>■ Pi, P 2 form a weak dual pair. Let h G C 1 (Pi). Let x(t) G S 
be a solution to Hamilton’s equations with respect to the Hamiltonian H = ho J 1 . Then 
J\ (x(t)) G Pi is a solution to Hamilton’s equations on P\ with respect to h, and J 2 (x(t )) is 
constant in time. 


Proof. Use proposition A.7 to show that fi(t) = Ji(x(t)) is a solution to Hamilton’s equations 
with respect to h. To verify that J 2 (x(i)) is constant, let v G kernel(Ji) be a vector over 
x(0) G S. This means that v is tangent to the level set of J\ at x(0). Moreover, H = h o J 1 
is constant on such level sets. Thus we observe 


0 = (dH(x(Q)),v) = lu (x, v ). 

Since v was an arbitrary element of kernel( J\) over x(0) we see that x G kernel( J\) u . Since 
Ji and Jo form a weak dual pair, this implies x G kernel(J 2 ). Thus we have found 

^J 2 (x(t)) =TJ 2 ■ x(t) = 0. 

□ 


Appendix B. Diagrammatic overview 


We present here a diagrammatic representation of some of the spaces used in the present 
paper. We begin by recalling a number of general results that hold for finite-dimensional Lie 
groups, before we indicate their relevance to the developments in the main text. 



Q 



q = H(q) = (p g (q 0 ) 
-•- 




Figure 6. The G action on Q induces a projection n. The group K also 
acts on G, and we assume that the two actions commute. 


• Let a Lie group G act on a manifold Q by the action </> : G x Q -A- Q, which we 
also write as = <f>(g, •)• If we fix a particular value q 0 G Q, we can construct 
a mapping n : G —> Q given by n(r/) = f> g (qo). Let us assume that the action is 
transitive, so that n is surjective. We denote by G qo the isotropy subgroup leaving 
q 0 invariant, that is, 

Gqo := {9 £ G\4> g (q 0 ) = go}- 

Note that n _1 (g) = gG qo for any g G G such that n(g) = g, and hence we can 
identify Q with G/G qo . Suppose a further Lie group, K, also acts on Q with group 
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action 0 : Q x K —>■ Q, which commutes with 0. This situation arises naturally, for 
instance, when K is a subgroup of G and, in turn, G qo is a normal subgroup of K. 
In that case, one can define the action 0 as 

(28) ip k (q) := Tl(gk) = 0 sfc (g o ), 

and check that 0 and 0 indeed commute: 

0 S ('0fc(g)) = 0s(</Vc(<?o)) = <i>sgk(qo) = 0fc(0s (<?))■ 

We refer to Figure [6] for a representation of the relevant spaces and maps. 

• The actions 0 and 0 on Q can be lifted to actions $ and T on the cotangent bundle 
T*Q in the usual manner (see Figure [7]). These cotangent lifted actions induce 
equivariant momentum maps J\ : T*Q —> g* and J 2 : T*Q —> £*, where g* and t* 
are the duals of the Lie algebras of G and K . Due to their equivariance, J\ and J 2 
are Poisson maps (where the duals of the Lie algebras are equipped with appropriate 
Lie-Poisson brackets). Since the actions 0 and 0 commute, the action T leaves level 
sets of ,/] invariant, and vice versa. This implies that J\ and J 2 are a weak dual pair, 
and if moreover T is transitive on the level sets of J\ and vice versa, then J\ and J 2 
are a proper dual pair, see [ GBVl 2 ! . Corollary 2.6]. 

• Let H : T*G Aftbea right-invariant Hamiltonian. This means that H is invariant 
with respect to the cotangent lift TR* of the multiplication from the right of G 
by itself. In particular, the reduced Hamiltonian h : g* —> M satisfies "H(a: g ) = 
hoTR*(a g ) for any a g G T*G , and the reduced dynamics in g* are of Lie-Poisson type. 
The momentum map J\ can be used to induce the so-called collective Hamiltonian 
H = ho J 1 on T*Q. Note that T*Q is a symplectic manifold, and that the symplectic 
(canonical) dynamics with respect to the collective Hamiltonian are mapped by (the 
Poisson map) J\ to the reduced dynamics on g*. Moreover, J 2 is conserved under 
the dynamics on T*Q. The conservation law follows from Noether’s theorem because 
Ji, and hence H, are left invariant by T (see Theorem A.9). Note that the elements 
of T*Q play the role of symplectic variables (or Clebsch variables in the sense of 
|MW83| L 

• The appeal of Clebsch variables is their symplectic nature. The symmetry of H 
with respect to T implies that reduced dynamics on T*G/H can be constructed 
by symplectic reduction. Note however that the resulting quotient manifold is not 
symplectic in general. 

• Note also that there is a symplectic diffeomorphism between T*Q and J _1 (0)/G, 


90) 


where J here is the momentum map associated with the cotangent lift of the action 
(from the right) of G qo on G , see [MMQ + 07[ Theorem 2.2.2], 

In translating the above facts to the case of interest in the present paper, one encounters 
technical subtleties to do with the infinite-dimensionality of SDiff(M n ). Nevertheless, bullet- 
by-bullct parallels can be recognized between the developments in the main text of the 
paper and the general results above, as we will discuss now. For simplicity, we restrict 
ourselves in what follows to the case of a single particle, the extension to N particles being 
straightforward. 

• Let G = SDiff(R n ), and let Q = Q[ k) be the space of single-particle fc-jetlets. We 
fix the point q 0 = (z, 1 , 0,... , 0) G Q[ k) corresponding to the Taylor expansion of the 
identity in SDiff(M n ) evaluated at z G (cf. Section 4.4). We let the projection n 
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T*G 



tr G 

K\ 


3 


X Gqo 

> e 




FIGURE 7. The cotangent lifted actions induce momentum maps J\ and J r 2 . 
The Hamiltonian "H, assumed to be right-invariant, is also shown, along with 
the reduced Hamiltonian H and the collective Hamiltonian H — h o J,. 


be given by Jet^, and hence the left action of an element <p G SDiff(M n ) on q G Q 
is 4> v (q) = Jeto p), for any p such that Jet^(p) = q. The role of the isotropy 
subgroup G qo is played by iso^. Let K = iso(z), and note that iso^^z) is a normal 
subgroup of iso(z) ([ DJR.13 , Proposition 4.1]). Hence, we can define a right action of 


iso(z) on given by (28), namely 

I’piq) = Jet f } (<p°p), 


for any <p such that Jet^(<p) = q. 

The cotangent lift of the left and right actions on Q lead to the momentum maps 

(k) (k) 

J r and Jp , explicitly computed in Section 4.4 and shown to be a weak dual pair 


in Proposition |4. 7 


associated with the collective Hamiltonian = h P)(y o J L 
tJ) 


(fc) 


In Theorem |4.8| we showed that the canonical Hamiltonian equations on T*Q 

lead to trajectories that 
are mapped, by J} K; , to solutions of Hamilton’s equations on 3C div (M n )*. Moreover, 
we showed that is a constant of motion. 

We briefly visit further symplectic reduction in Appendix [Ej For more on this topic 
we refer t o [DJR13 . CHJM14 J. 

In Section 


4.5 


(18), namely, 


we constructed the mapping i : T*Q ^ -A T* SDiff(M n )/iso^ fc ^ given by 


i(q,p) = [TR;^jP(q,p)]. 

As we briefly mentioned towards the end of that section, if we denote by J c 
T*SDiff(M n ) —y iso^J(z)* the momentum map associated to the cotangent lift of 
the right action of iso^z) on SDiff(M n ), then one can show that i(T*Q = 


(k) 


(*) N-l 


( jw ) 

V^conv ) 


(0)/iso^(z), as suggested by |MMQ + 07l. Theorem 2.2.2], 
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Appendix C. Multi-indices 

A multiset is a set with some notion of multiplicity [Bli89] . In this paper, a multi-index 
on M n is a multiset of elements derived from the generating set {1,... ,n}. Heuristically, a 
multi-index is just a “bag of marbles” each of which comes in n “colors”. Given two multi¬ 
indices a G bag* (n) and f3 G bag J (n) one can create the multiset union a U /3 G bag 1+3 (n) 
by collecting the marbles of a and f3 into a single bag. Given integers bi,... ,bj G M”, we 
can define the unique multi-index f3 = [bi,... ,bj] G bag 3 {n) obtained by collecting bi,... ,bj 
into a bag. Given these conventions, we can denote the partial differential operator db 1 ...b n 
by dp. Moreover, the notion of equivalence of mixed partials is expressed by the equivalence 
d a dp = d aU p = dpd a . The cardinality of the multi-index j3 is denoted \f}\ and is given by 
the number of marbles in the bag. Thus the order of the partial differential operator dp is 
|/31. We denote the space of fc-th order partitions of a multi-index by II(a, k). Rather than 
defining all this formally, we will compute an example and refer to [Jacl4| for the formal 
definitions. 

We can consider the integers 1 2 and 1, and the partial differential operator d\ 2 i- The 
associated multi-index is just [1,2,1]. This multi-index is equivalent to the multi-index 
[1,1,2] and [2,1,1], We say that it contains the elements 1 and 2. Because it contains 
‘1’ two times, we say that the multiplicity of 1 is 2. The multiset of 2-fold partitions is 
n([1, 2,1], 2), and consists of three multiset-partitions 

[[1],[2,1 ]],[[1,1],[2]],[[1 ,2 ],[1]]. 

Note that the first and the third partition correspond to the same multiset. The cardinality 
of II([1, 2,1], 2) is 3, although it only has two distinct elements (one with a multiplicity of 1, 
and another with a multiplicity of 2). 

Appendix D. The dual space to divergence free vector fields 

In this section we will provide a terse and incomplete characterization of the dual space 
of divergence free vector fields. First let us characterize the dual space of the space of all 
vector fields (with “proper” decay). Let £(M n ) be the space of vector fields which decay 
at infinity in such a way that X(l n ) is Frechet. Viewing X(M n ) as a subspace of functions 
from M n to M n we can view its dual as a space of distributions. That is to say, given any 
m G 3C(M n )* we may write m as a tensor product p ® p where p is a distribution (perhaps a 
measure) on M n and p is a covector field (i.e. a one-form). Conversely, given any p G fT(M n ) 
and distribution p G V(R n ) we may form the tensor product p ® p. The object p ® p is 
identified as an element of X(M n )* through the pairing 

(p®p,u) := / (p(x),u{x))p. 

J R" 

where (p(x),u(x)) is the function on M n obtained by pairing the covector p(x) G with 

the vector u(x) G T x M n . If we restrict ourselves to the case of divergence free vector fields, 
we need to quotient the dual space appropriately. In particular, we see that the annihilator 
of 3Cdiv(R n ) as a subspace of 3C(M n )* is 

(X div (M n ))° := {m G £(M n )* | (m,u) = 0, Vw G 3C div (M n )} 

:= closure(jj <8) da; G 3C(M n )* | p is a closed one-form}. 
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where cLr is the canonical volume form on M” and we have used the fact that the gradient 
fields and the harmonic vector fields are L 2 -orthogonal to the divergence free vector fields. 
The dual space £div(®^ n )* is identical to the quotient space X(M n )*/(Xdi V (K n ))°- hi other 
words, we may view a me £div(B^ n )* as an object of the form p ® p modulo Xdi V (K n )°- 
In the text we will typically not mention £ di v(M n )° explicitly, and simply identify m with 
p <E> P- This is a harmless identification as long as we do not pair it with a non-divergence 
free vector field. 


Appendix E. Equations of motion for 1-jetlets 


The equations of motion are expressible as Hamiltonian equations on T*Q$ in canonical 
variable (g(°),p(°), gl 1 )^! 1 )). However, it is more efficient to express the equations of motion 

in the non-canonical variables (q,p,p) where q a = qll\ p a = pi°' > and p a = [q^YilP^]/ f° r 
a — 1,..., N. The Hamiltonian in these coordinates is 

H(q,p, p) = ip„(Py K v (q* - qb ) - Pai[Pb]j k 9kK tJ (q a - q b ) 

-l^WbV^^-qb), 

where K^(x) = 6^ e - ^^ 2 / 2cr “. Hamilton’s equations are then given in short by 


(29) 

(30) 

(31) 

(32) 


q = 

p = 


dH 

dp 

dH 

dq 


£ 

/j 


dH 

dp 

— adg(p) 


where ad* refers to the coadjoint operator on SL(n). More explicitly, equation (29) is given 
by 

= PbjK lJ (q a - q b ) - [p^ ] ] k d k K i: >{q a - q b ) 


equation (30) is given by the sum 


• _ rjlOO , rji 01 I rpll 

fat - L ai ' - L ai ' - 1 ai i 


where we define the three terms in this sum as 

T™ = -PakPbAK kj (q a ~ q b ) 

T al = (Pai[p>b ) ] j k - Pbi[p { a ] }j k )dkiK h (q a - q b ) 

Til = AAJl^\ k dik*K m \q a - q b ). 


Next, we calculate 

which allows us to 


the quantities £ = dH/dp for k — 1, 2 of equation (31) to be 
[£a]\ = PbkdjK lk (q a - q b ) - [pAdjkK^qa - q b ), 


compute p in equation (32) as 

W/ = MlePh - kUBl’A 
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The dynamics in terms of the original variables (q^°\q^) with q^ e SL(n) are obtained 

by integrating the reconstruction equations [q^Yj = [£]*fc(V 
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